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Abstract. The so-called internal shock model aims to explain the light-curves and spectra pro- 
duced by non-thermal processes originated in the flow of blazars and gamma-ray bursts. A long 
standing question is whether the tenuous collisionless shocks, driven inside a relativistic flow, are 
efficient enough to explain the amount of energy observed as compared with the expected kinetic 
power of the outflow. In this work we study the dynamic efficiency of conversion of kinetic-to- 
thermal/magnetic energy of internal shocks in relativistic magnetized outflows. We find that the 
collision between shells with a non-zero relative velocity can yield either two oppositely moving 
shocks (in the frame where the contact surface is at rest), or a reverse shock and a forward rar- 
efaction. For moderately magnetized shocks (magnetization a ~ 0.1), the dynamic efficiency in a 
single two-shell interaction can be as large as 40%. Hence, the dynamic efficiency of moderately 
magnetized shocks is larger than in the corresponding unmagnetized two-shell interaction. We find 
that the efficiency is only weakly dependent on the Lorentz factor of the shells and, thus internal 
shocks in the magnetized flow of blazars and gamma-ray bursts are approximately equally efficient. 
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1. INTRODUCTION 

Internal shocks (ISs) [1] are invoked to explain the variability of blazars [e.g., 2, 3] 
and the light curves of the prompt phase of gamma-ray bursts (GRBs) [4, 5, 6]. A long 
standing concern in this model is the question whether it is efficient enough to explain 
the relation between the observed energies both in the prompt GRB phase and in the 
afterglow [see, e.g., 7, 8, 9, 10]. To assess the efficiency of the internal shock model, 
most of the previous works focus on the comparison between the observed light curves 
and the model predictions employing a simple inelastic collision of two-point masses 
[7, 11, 12, 13, 14]. Less attention has been paid to the hydrodynamic effects during the 
shell collision [but see, 9, 15, 16, 3, 17]. 

The ejecta in GRBs and blazars may be significantly magnetized, particularly if 
they originate from a Poynting-flux-dominated flow [18]. Forming shocks in highly 
magnetized media is challenging since the Alfven speed approaches the speed of light. 
Therefore, to account for the observed phenomenology it is necessary to address how 
efficient can be shock dissipation of the internal collisions in arbitrarily magnetized 
flows. This question has been partly considered by a few recent works [e.g., 19, 20], 
and only very recently has been studied extensively [21]. 

The base for a study the efficiency of internal collisions is the determination of the dy- 
namic efficiency of a single binary collision, i.e., the efficiency of converting the kinetic 



energy of the colliding fluid into thermal and/or magnetic energy. Thus, we model ISs 
as results of the collision of (magneto-)hydrodynamic shells of plasma with a non-zero 
relative velocity. The contact surface, where the interaction between the shells occurs, 
can break up either into two oppositely moving shocks (in the frame where the contact 
surface is at rest), or into a reverse shock and a forward rarefaction. The determination 
of whether one or the other possibility occurs is computed by estimating the invariant 
relative velocity between the fastest and the slowest shell, i.e., by solving the relativis- 
tic magnetohydrodynamic (RMHD) Riemann problem posed by the piecewise uniform 
states given by the physical quantities on the two interacting shells (Sect. 2). In Section 3 
we define precisely the notion of dynamic efficiency, both for shocks and rarefactions. 
We perform a parametric study of the binary shell collision dynamic efficiency in Sec- 
tion 4. Finally, the discussion and conclusions are listed in Section 5. 



The interaction between parts of the outflow with varying properties can be modeled in 
terms of Riemann problems, i.e. relativistic magnetohydrodynamic initial-value prob- 
lems with two constant states separated by a discontinuity. This approach allows us to 
use our models to sample very finely a large parameter space. The approach is also jus- 
tified by the fact that the flow is cold and ultrarelativistic, so that its lateral expansion is 
negligible. Thus, a description of the interactions assuming planar symmetry suffices to 
compute the dynamic efficiency of such interactions. 

In the following we use subscripts L and R to denote properties of the (faster) left and 
(slower) right state, respectively. We normalize the rest-mass density p to pR, the energy 
density to pRC 2 (c is the speed of light) and the magnetic field strength to cy/4Wp~R. 

For the initial thermal pressure of both states we assume that it is a small fraction 
of the density, pi = xpL and pr = %, and that magnetic fields are perpendicular to the 
direction of the flow propagation. The remaining parameters determining the RMHD 
Riemann problem are: the density contrast pi, the Lorentz factor of the right state Tr, 
the relative Lorentz factor difference Ag := (Ti — Tr) /Tr, and the magnetizations of left 
and right states, Ol '■= B\/ (r|(l + Ag) 2 pi) and Or := B 2 R /T 2 R , where Bl and Br are the 
lab frame magnetic field strengths of left and right states, respectively. Furthermore, we 
define the total (thermal + magnetic) pressure, the total specific enthalpy and the total 
energy density, respectively, as 



where £ denotes the specific internal energy. 

The typical structure of the flow after the break-up of the initial discontinuity consists 
of the two initial states, and two intermediate states separated by a contact discontinuity 
(CD) -see [22]. The total pressure and velocity are the same on both sides of the CD. 
The quantity o/p is uniform everywhere, except across the CD, where it can have a 



2. RMHD RIEMANN PROBLEM 




(1) 



(2) 



e* :=p(l + 8) + ^, 



(3) 



jump. We denote the total pressure of intermediate states p* s , and rest-mass density left 
and right of the CD as p S L and p S R l . 

One of the key steps in solving a Riemann problem is to determine under which 
conditions shocks can form. This happens when the Lorentz-invariant relative velocity 
between the left and right states, measured in the frame of the CD, v/ r := (v/ — v r )/(l — 
v/v r ) is larger than the limiting value [21] 



Generally, the quantity (v/ r )2s can be only determined numerically. If (v/ r ) < (v/ r )2s, 
a single shock and a rarefaction emerge from the initial discontinuity. 



To study ISs we idealize interactions of parts of the outflow moving with different 
velocities as collisions of homogeneous shells. In our model the faster (left) shell catches 
up with the slower (right) one yielding, in some cases, a pair of shocks propagating in 
opposite directions (as seen from the CD frame). In order to cover a wide range of 
possible flow Lorentz factors and shell magnetizations, we assume that initially, the flux 
of energy in the lab frame is uniform and the same in both shells (see [21]). We then 
compute the break up of the initial discontinuity between both shells using the exact 
Riemann solver developed by Romero et al. [22], assuming an ideal gas equation of 
state (EoS) with an adiabatic index 7= 4/3. 

Efficiency of energy dissipation by a shock. To model the dynamic efficiency of 
energy dissipation we follow the approach described in [21]. By using the exact solver 
we determine the existence of shocks and (in case one or two shocks exist) obtain the 
hydrodynamic state of the shocked fluid. We use subscripts 5,L and S,R to denote 
shocked portions of left and right shells, respectively. In the following we treat the 
efficiency of each shock separately. 

Reverse shock. To compute the efficiency we need to compare the energy content of 
the initial (unshocked) faster shell with that of the shocked shell at the moment when RS 
has crossed the initial shell. Assuming an initial shell width Ax, we define total initial 



In the context of ISs, if the flow is ultrarelativistic in the direction of propagation, the velocity compo- 
nents perpendicular to the flow propagation should be negligibly small and, hence, they are set up to zero 
in our model. If such velocities were significant, appreciable changes in the Riemann structure may result 
as pointed out in Aloy and Rezzolla [23] or Aloy and Mimica [24]. 




(4) 



3. ENERGY DISSIPATION EFFICIENCY OF ISS 



kinetic, thermal and magnetic energy [21] 



E K (r,p,Ax) := r(r-l)pAx, 
E T (r,p,p,Ax) := [(pe + p)T 2 -p]Ax, (5) 
E M (T,p,a,Ax) := (r 2 -1/2) pa Ax. 

When the RS crosses the whole initial shell, the length of the compressed shell (i.e., the 
fluid between the RS and the CD) is ClAx, where £l := (vcd — vs,l) I ( v l — vs.l) < 1> 
and vcd an d vs,l are velocities (in the lab frame) of the contact discontinuity and the RS, 
both obtained from the solver. We normalize the energies taking Ax = 1, and define the 
dynamic thermal and magnetic efficiencies, i.e., the fraction of the initial energy that the 
RS has converted into thermal and magnetic energy, respectively, as 

ET(^S,L,Ps,L,PS,L,CL)-E T (TR(l+Ag),p L ,XpLA) ,„ 

£t,l := , (o) 

Eo 

Em(^s,l,Ps,l, Os,l, Cl) -E M (T R (l+Ag),p L , a L , 1) 

£m,l := , (/) 

Eo 

where Eo is the total initial energy of both shells 

E := E K (r R (l+Ag),p L , l) + E K (T R , 1, 1) + E T (r R (l+Ag),p L ,xp L , 1)+ 

(8) 

E T (T R ,l,x, l)+E M (T R (\ + Ag),p L ,o L , 1)+E M (T R , I, Or, 1) . 

Forward shock. In complete analogy, we define the thermal and magnetic efficien- 
cies for the forward shock, 

E T (r SjR ,ps, R ,Ps,R,£R)-E T (r R ,l,X,l) m , 

£t,r-= , (y) 

Eo 

Em(Fs,r,Ps,r,Gs.r, Cr)-E m (T r , 1,o r , 1) 
£m.r '■= ^ , (10) 

where £ R := (vs. R — vcd) / {vs,R — vr) < 1, and V5 ; ^ is the lab frame velocity of the FS. 

Combining equations (6), (7), (9) and (10) we define the dynamic thermal and mag- 
netic efficiency of ISs 

£t ■= £t,l + £t,r, £m-=£m,l + £m,r- (11) 

Efficiency of energy dissipation by a rarefaction. In a rarefaction there is a net 
conversion of magnetic and/or thermal energy into kinetic energy, thus the net dynamic 
efficiency produced by a rarefaction, defined as in e.g., Eq. (7), should be negative. 
Therefore, it is possible that in some cases the total (left plus right) thermal or magnetic 
efficiency (Eqs. (11)) becomes negative. However, this situation does not correctly 
model the fact that, in cases where a single shock exists, it is still able to radiate 
away part of the thermal or magnetic energy behind it. Thus, we set £t,l = £m,l = 
(£t,r = £m,r = 0) if the reverse (forward) shock is absent. 
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FIGURE 1. Panel (a): contours of the total dynamic efficiency Ej + Em (Eq. (1 1)) in the GRB regime 
(Yr = 100, Ag — 1; left panel) for different combinations of (ol,Or). Contours indicate the efficiency 
in percent and their levels are 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12 and 13. In the region of the parameter 
space above the dashed line there is no forward shock, while the reverse shock is always present for 
the considered parametrization. Colors: magnetic efficiency Em in percent. Panel (b): The color scale 
indicates the value of the maximum total dynamic efficiency (in percent) as a function of the parameter 
pair (Ag,As). The values of the rest of the parameters are fixed to Fr = 100, and % = 10~ 4 . Contours: 
magnetization of the slowest shell: Or = 0.1, 0.5, 1, 5 and 10. 



4. PARAMETRIC STUDY OF THE DYNAMIC EFFICIENCY 

Next we study the dynamic dissipation efficiency in the process of collision of cold, 
magnetized shells. The shells are assumed to be cold because in the standard fireball 
model, almost all the internal energy of the ejecta has been converted to kinetic energy 
before ISs happen. If the ejecta were accelerated by magnetic fields, then the flow is 
cold (i.e., % «C 1) all the way from the beginning to the internal shock phase. Thus, 
we fix x = 10 4 <C 1, to model initially cold shells. We set Ag = 1 as a reference 
value, and consider the dependence of the total dynamical efficiency (£j + £m) on the 
individual magnetization of each shell Ol and Or Fig. 1(a). For the purpose of this study 
we take Tr = 100, which is a typical value of GRB outflows. Smaller values, = 10, 
more representative of blazar outflows, do not yield significant differences in the total 
dynamical efficiency [21]. 

The maximum efficiency (where £j + £m > 0.13) is attained for (Or, o£) w (0.2, 1). 
In the region above the dashed line of Fig. 1(a) the FS is absent and, thus, since only 
the RS dissipates the initial energy, the efficiency slightly drops. However, the transition 
between the regime where the two shocks operate or only the RS exists is smooth. The 
reason being that the efficiency below the separatrix of the two regimes and close to it is 
dominated by the contribution of RS. As expected, when either Ol or Or approach low 
values, the dynamic efficiency ceases to depend on them. 

We illustrate the differences in terms of dynamical efficiency by choosing three 



characteristic models in the parameter space (letters A, B and C; Fig. la). Model A 
corresponds to a prototype of interaction between non-magnetized shells (oi = Or — 
10~ 6 ), model B is located at the maximum efficiency (ol = 0.8, Or = 0.2), and model 
C corresponds to the case when the FS is absent (ol = 1, Or = 10 2 ). All three models 
have a substantial dynamic efficiency, but there is a qualitative difference among them. 
In model A, ISs dissipate kinetic to thermal energy only (thermal efficiency). In model 
B shocks mainly compress the magnetic field (magnetic efficiency) and dissipate only a 
minor fraction of the initial kinetic and magnetic energy to thermal energy. Finally, in 
the model C only the RS is active. 

Values Ag < 2 between adjacent parts of the flow are motivated by the results of 
numerical simulations of relativistic outflows [e.g., 26, 27, 28, 29, 30, 31, 32, 33] (but 
see, e.g., [34], who find Ag ~ 1 — 19 appropriate to model Mrk 421). This adjacent flow 
regions can be approximated as pairs of shells whose binary collision we are considering 
here. However, it has been confirmed by several independent works (e.g., [7, 8, 9, 15]) 
that, in order to achieve a high efficiency (more than a few percent) in internal collisions 
of unmagnetized shells, the ratio between the maximum (r max ) and the minimum (r m i n ) 
Lorentz factor of the distribution of initial shells should be r max /r m i n > 10. 

In view of these results, we have also made an extensive analysis of the dependence 
of the dynamic efficiency on Ag. Since we are also interested in evaluating the influence 
of the magnetic fields on the results, we define a new variable As = ( 1 + Ol) /( 1 + Or), 
and plot (Fig. lb) the value of the maximum efficiency reached for every combination 
(Ag,As) and fixed values of the rest of the parameters (T R = 100, x = 10~ 4 ). We find 
that the maximum total dynamic efficiency grows (non-monotonically) with increasing 
Ag (Fig. lb), in agreement with the above mentioned works for unmagnetized shell 
collisions. Indeed, a large value Ag J> 10 yields dynamic efficiency values ~ 40% if both 
shells are moderately magnetized (Or ~ Ol ^0. 1). Nevertheless, the amount of increase 
of efficiency with Ag depends strongly on As. For |As| > 1, corresponding to cases where 
the slower shell is highly magnetized (Or > 4), the maximum dynamic efficiency is 
almost independent of Ag; while for | A$| <, 1, the maximum dynamic efficiency displays 
a strong, non-monotonic dependence on Ag. 

It is remarkable that values 5 <, As <^ 100 yield dynamic efficiencies in excess of 
~ 20%, regardless of the relative Lorentz factor between the two shells. In this region of 
the parameter space the maximal dynamic efficiency happens when Or > 10, Ol > 50, 
and the total dynamic magnetic efficiency dominates the total dynamic efficiency. 

5. DISCUSSION 

In this work we estimate the dynamic efficiency of conversion of kinetic-to- 
thermal/magnetic energy in collisions (ISs) of magnetized shells in relativistic outflows. 
A fundamental difference between the internal collisions in magnetized and unmagne- 
tized outflows is the fact that in the former case not only shocks but also rarefactions can 
form. Thus, one would naturally expect a reduced dynamic efficiency in the magnetized 
case. However, we find that the dynamic efficiency reaches values ~ 10% —40%, in 
a wide range of the parameter space typical for relativistic outflows of astrophysical 
interest (blazars and GRBs). Thus, the dynamic efficiency of moderately magnetized 



shell interactions is larger than in the corresponding unmagnetized case. This is because 
when the shells are moderately magnetized, most of the initial shell kinetic energy is 
converted to magnetic energy, rather than to thermal energy. 

The difference in efficiency between flows in blazars (Tr = 10) and GRBs (Tr = 100) 
is very small if Ag is fixed. From theoretical and numerical grounds, values of Ag ~ 1 
seem reasonable, and Ag = 1 has been taken as a typical value for both blazars and 
GRB jets, which brings maximum efficiency when the magnetizations of the colliding 
shells are (<Jl,Or) ~ (1,0.2). Larger dynamic efficiency values ~ 40% are reached if 
4? ^ 10 an d | As | <^0.1, corresponding to cases where the magnetization of both shells 
is moderate (Gr ~ Ol < 0.1). 

In the limit of low magnetization of both shells, the kinetic energy is mostly converted 
into thermal energy, where the increased magnetic energy in the shocked plasma is only 
a minor contribution to the total dynamic efficiency, i.e., £j <C £m- Here we find that 
as the magnetization of the shells grows, the roles of £? and £ M are exchanged, so 
that £j < £m (at the maximum dynamic efficiency £j ~ 0.1 £m)- If the magnetization 
of both shells is large, the dynamic efficiency decreases again because producing shocks 
in highly magnetized media is very difficult. All these conclusions are independent on 
the EoS used to model the plasma [21]. 

The comparison of our results with previous analytic or semi-analytic works is not 
straightforward, since, generally, they do not compute the (magneto-)hydrodynamic 
shell evolution, and they include multiple interactions of a number of dense shells. The 
bottom line in these previous works is that internal collisions of unmagnetized shells 
can be extremely efficient (~ 100%; [8]) if the spread of the Lorentz factor (i.e., the 
ratio between the Lorentz factor of the faster, r max , and of the slower r m i n shell in 
the sample) is large (r max /T m j n = 10 3 ; e.g., [7, 35]). A more moderate spread of the 
Lorentz factor T max /T m i n = 10, yields efficiencies ~ 20%. These high efficiencies result 
after a large number of binary collisions. Since more than a single collision is included, 
the kinetic energy remaining in the flow after the first generation of collisions, can be 
further converted into internal energy as subsequent generations of collisions take place. 
In contrast, we find that moderate magnetizations of both shells (o ^0.1) and Ag J> 10 
(which would roughly correspond to r max /r m i n = 9) are enough for a single binary 
collision to reach a total dynamic efficiency of ~ 40%. 

We point out that the energy radiated in the collision of magnetized shells is only a 
fraction, f r < 1 [e.g., 36, 37, 8] of the energy dynamically converted into thermal or 
magnetic energy.Hence, a single binary collision between moderately magnetized shells 
may yield radiative efficiencies ~ 0.4/ r . Therefore, (1) binary collisions in relativistic 
magnetized flows are efficient enough, from the dynamical point of view, to be a valid 
mechanism to dissipate the bulk kinetic energy of relativistic ejecta, and (2) the main 
restriction on the radiative efficiency comes from the radiation mechanism setting the 
limiting factor f r . The estimated dynamic efficiency in the binary collision of magne- 
tized shells will be considered in a future work by computing the radiative efficiency 
using the method devised in [38]. 
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